The Two-Impurity Anderson Model at Quantum Criticality 
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We propose a realization of the two-impurity Anderson model in a double quantum-dot device. 
When charge transfer between the dots is suppressed the system exhibits a quantum phase transi- 
tion, which is controlled by a surface of non-Fermi liquid fixed points parameterized by the charge 
valences of the dots. Employing conformal field theory techniques, we identify the scaling exponents 
that govern transport and thermodynamics close to criticality. We also determine the dynamical 
exponents that set the time scale for the buildup of the non-Fermi liquid state after the system is 
suddenly shifted into the critical region, e.g. by a change of a nearby gate voltage. 
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I. INTRODUCTION 

The physics of two interacting quantum impurities em- 
bedded in an electron gas is an important and much stud- 
ied problem. It has bearings on the theory of quantum 
critical phenomena [l[ and turns up in a host of other 
problems, including the study of heavy fermion physics 
[2| , cluster implementations of dynamical mean- field the- 
ory [|[ , and the quest for a solid-state quantum computer 

!■ 

When the impurities carry spin-1/2 magnetic mo- 
ments, the problem becomes that of the two-impurity 
Kondo model ("TIKM") 0]. This model supports two 
competing energy scales: an Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interaction ~ K driven by the conduc- 
tion electrons — which, for K > 0, tends to lock the 
magnetic moments into a singlet — and a Kondo tem- 
perature Tk below which the electrons try to screen 
the moments. When the two scales become compara- 
ble, there is a crossover between the RKKY-singlet and 
Kondo screened phases @ . If the model is endowed with 
a special kind of electron-hole symmetry, the crossover 
sharpens into a second-order phase transition, controlled 
by a non- Fermi liquid fixed point Q. However, since 
the required symmetry is not expected to be found in 
any known material, the experimental observation of the 
quantum critical state has been seen as rather unlikely. 

In a recent work, Zarand et al. [1] suggested that an ex- 
perimentally controlled approach to the TIKM quantum 
critical state may in fact be achieved by using a special 
type of double-quantum dot device, realized in a gated 
semiconductor heterostructure. With a design where the 
two dots are connected to two separate leads, and RKKY- 
couplcd via a magnetic insulator, the quantum critical 
state was predicted to be insensitive to electron-hole sym- 
metry breaking. 

In this paper we address the question what happens if 
the proposed device is operated in a regime where charge 
is allowed to fluctuate on the dots (but with charge trans- 
fer between the dots still being suppressed)? An answer 
to this question is important for validating future ex- 
perimental tests of the proposal in Ref. [g], as well as 



for characterizing quantum criticality in the presence of 
a local perturbation, effectively generated by processes 
at a higher energy scale. To this end we study the two- 
impurity Anderson model ("TIAM") 0], with the two im- 
purities connected to two separate leads. Applying an ex- 
tension of the boundary conformal field theory approach 
(BCFT) that has been used for the TIKM @], we are able 
to explore the quantum critical properties of this model 
non-perturbatively. We find that the charge fluctuations 
generate a surface of unstable fixed points, connected to 
that of the TIKM fixed point by two marginal opera- 
tors, but with thermodynamics and transport properties 
affected only in system-specific amplitudes. In contrast, 
the dynamical exponents that govern the approach to 
the non-Fermi liquid state after the system has suddenly 
been shifted into the critical region sensitively depend on 
how much charge is localized on each dot. 

The setup we consider is depicted in Fig. [TJ The two 
infinite leads are coupled to their respective dots via tun- 
nel junctions, allowing the electrons to hop from [into] 
lead i into [from] dot i with amplitude Vi (i = 1,2). The 
finite auxiliary electron reservoir (a large quantum dot 
with a quasi-continuous density of states p) is coupled to 
both dots also via tunnel junctions, with tunneling rates 
~ V A a. 




Lead 1 



Lead 2 



FIG. 1: (Color online) The physical system that we study in 
this article. The different V are tunneling rates. The auxiliary 
reservoir is operated in the Coulomb blockade regime, where 
charge transfer between the reservoir and dots (and thus also 
between the two leads) is strongly suppressed. 

First consider the case V\ = V<x = 0. If we tune nearby 
gate voltages such that each dot has spin 1 /2, with charge 
transfer between the dots and the auxiliary reservoir sup- 
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pressed by Coulomb blockade [10( , the setup is described 
by the TIKM, with Hamiltonian 

#tikm = H kin + ^2 Ja.iVi ■ (?a{Ri) + Kai ■ er 2 , (1) 

i=l,2 

where J as ~ Vaa/Ua (i = 1,2), and Ua being the 
Coulomb blockade energy of the reservoir. Here = 
d\ a T a pdip and cr A (x) = c^x^a^c^x), where d\ a cre- 
ates an electron with spin a on dot i, c' a creates an elec- 
tron with spin a in the auxiliary reservoir, and Hi and 
i?2 denote the locations of the two dots. T a p are elements 
of the vector of Pauli matrices, and all repeated spin in- 
dices are summed over, i/kin is the kinetic energy of the 
electrons in the auxiliary reservoir. The RKKY interac- 
tion with coupling K ~ Ja,iJa,2 cosip/D is generated in 
second order in the Kondo couplings J as, with D the 
band width of the reservoir and where ip ~ R\ — i?2 con- 
trols the sign of the interaction. When K » Ta, with 
Ta being the Kondo temperature associated with the 
auxiliary reservoir, the antifcrromagnctic RKKY inter- 
action Kcr\ ■ (T2 dominates the direct Kondo interactions 
0, HH , and we are left with only the first and last term in 
the Hamiltonian in Eq. ([1]). Here Ta = min(TA i, Ta,2), 
with Tas ~ D exp(— \j1pJ A ,i)- If we now turn on V1.2, 
charge may fluctuate between each of the dots and their 
respective leads, but not between the two dots. By this 
construction the system can reach a quantum critical 
phase, protected by the suppression of charge transfer 
between the two leads @. As we shall see, this phase 
is not fragile against anisotropics in the couplings Vi, as 
long as charge transfer between the leads is suppressed 
(which is achieved via Coulomb blockade, as described 
above). This is an important feature of the proposed 
device, as fine-tuning of the couplings is not necessary 
for accessing the critical state. With the condition that 
K » Ta, the Hamiltonian modeling the two dots and 
the electrons in the leads is that of the two-impurity An- 
derson model, 

-Htiam = Kcr 1 • <T2 + Hi + H%, (2) 

where 

H = i/ k jn,i + (Vic\ a {Ri)d ia + /i.e.) 

+ Uiid^di^d^dn) + €id\ a di a , t = 1,2. (3) 

Here H\^ n i is the kinetic energy for the electrons in lead 
i and c\ a creates an electron with spin a in lead i. U{ 
is the corresponding charging energy and is the chem- 
ical potential that is measured from the common Fermi 
level of the leads. If |e.;| is of the same order of magni- 
tude as the level width I\ ~ V.?, its singly occupied sub- 
spaces are hybridized with the doubly and unoccupied 
subspaces ("mixed valence regime"). At the electron- 
hole symmetric point ej = —Ui/2 (center of the Coulomb 
blockade valley of dot i) the net charge on the impurity 
is m = {d\ a d ia ) = 1, while n j; < 1 for larger values of tj. 
The rij > 1 case need not be considered separately as it 
is related to rii < 1 via an electron-hole transformation. 



II. IDENTIFYING THE CRITICAL THEORY 

As was first realized by Affleck and Ludwig [l2| . the 
critical properties intrinsic to a quantum impurity prob- 
lem are encoded by a boundary conformal field theory 
(BCFT). More precisely, by representing the impurity- 
electron interaction at criticality by a conformally invari- 
ant boundary condition, the electrons can be described 
by a BCFT. The method is well-documented in the lit- 
erature [l4|, and here we only sketch the basics so as to 
lay down notation and conventions. 

Any application of BCFT requires that the electrons 
are effectively confined to a one-dimensional (ID) geom- 
etry. In the present case, the point-like interaction of 
Eq. © affects only the I = components of the angular- 
momentum decomposition of the electron fields in each 
lead. The system can thus be described by a trivial non- 
interacting two-dimensional or three-dimensional theory 
(involving all higher angular-momentum states) and an 
interacting ID theory (involving only the I = states), 
where the interesting physics happens and on which we 
focus exclusively. With this, the first step in the BCFT 
approach is to set up a framework where the impurity- 
electron interaction can conveniently be traded for a 
boundary condition. To this end one may exploit a well- 
known result from conformal field theory that connects 
the energy levels of the BCFT on a finite strip (with com- 
plex coordinates {w = u+iv | —00 < u < 00, ^ v ^ £}) 
to the (boundary) scaling dimensions of operators in the 
upper complex half-plane {z = r + ix \ x ^ 0}. (Here r is 
the Euclidean time and x is the spatial coordinate of the 
ID model, with x — being the location of the impuri- 
ties (l5|.) By imposing a conformally invariant boundary 
condition, call it A, at the edges v = and v = £ of the 
strip, and then mapping it onto the semi-infinite plane us- 
ing the conformal transformation z = exp(irw / £) (imply- 
ing boundary condition A at x = 0), one obtains a con- 
nection between the finite-size energy spectrum on the 
strip and the boundary scaling dimensions in the semi- 
infinite plane. With Eq as the ground-state energy, one 
has 



where {E n } is the spectrum of excited energy levels in 
^ v ^ £, and {A„} is the spectrum of boundary scal- 
ing dimensions. The boundary scaling dimensions A„ 
control the asymptotic behavior of the autocorrelation 
functions close to the boundary, 

(lfi n (T,x)lfi n (0,x)) - (cp n (T,x))((p n (0,x)) ~ r - 2A " 

\t\»\x\, (5) 

where ip n is a conformal field with the property that the 
state <p n \0) has energy E n , with |0) being the ground 
state on the strip. The autocorrelation functions in 
turn determine how the presence of the boundary ( alias 



the impurity-electron interaction) influences the low- 
temperature thermodynamics and transport in its neigh- 
borhood (treating Euclidean time as an inverse temper- 
ature). The problem has thus been reduced to deter- 
mining the finite-size spectrum, Eq. ((4]), of the the- 
ory. First, however, the specific boundary condition A 
that represents the electron-impurity interaction has to 
be nailed down. This is a nontrivial task, and we first re- 
view how to carry it out in the somewhat simpler case of 
the TIKM 0] (which corresponds to the integer-valence 
limit Ui = -2e t -> oo (i = 1, 2) of the TIAM). 

The two separate channels of ID electron fields (rep- 
resenting the conduction electrons in leads 1 and 2) each 
have a U(l)-charge and an SU(2)-spin symmetry. Using 
non-Abelian bosonization, the fermionic theory can then 
be rewritten exactly in terms of two U(l)-charge bosons 
and two SU(2)i spin bosons, where k = 1 is the level 
of the SU{2)k algebra under which the bosons transform 
fl3| . To incorporate the fact that the Kondo interaction 
couples the two spin bosons and, therefore, breaks the 
SU(2)i<8> SU(2)i symmetry down to the diagonal sub- 
group ( "total spin" ) , the spin bosons are expressed in 
terms of a single SU(2)2 spin boson and — as required by 
the coset construction [13[ — an Ising field. The states in 
any conformal field theory organize into conformal tow- 
ers, corresponding to the irreducible representations of 
the symmetry group(s) of the theory. In the case of the 
TIKM (as well as the TIAM) the conformal towers are 
labeled by the quantum numbers Qi,2 S Z (net charge 
in respective lead with respect to the filled Fermi sea), 
j e {0, 1/2, 1} (total SU(2) 2 spin), and A G {0, 1/16, 1/2} 
(corresponding to the Ising fields I, a, e) . The identifica- 
tion of spin conformal towers is shown in Table |TJ 



TABLE I: Identification of spin conformal towers from the 
coset construction. 

SU(2) 1 ® SU(2) 1 SU(2) 2 g Ising 



(0)i <g> (0)i 
(l/2)i® (0)i 
(l/2)i ® (l/2)i 



[(0) 2 ®I]©[(l) 2 ®e] 
(l/2) 2 ®a 

[(0)2 06] ® [(1)2® I] 



All interactions between the bosons (originating from 
the Kondo- and RKKY-interactions in the electronic the- 
ory) are absorbed into a conformally invariant bound- 
ary condition so that the boson dynamics is governed 
by a free Hamiltonian. The boundary condition is en- 
coded in a set of gluing conditions for combining the 
conformal towers of the effective (conformal field) the- 
ory. The gluing conditions in turn constrain the finite- 
size energy spectrum from which the dimensions of the 
scaling operators — and thus the critical scaling of ob- 
servables — are extracted. In the present case, it has 
been shown that the particular gluing conditions that 
describe the impurity-electron interaction at the critical 
point — which is known to occur for some intermediate 
value of the RKKY coupling, K ~ T K = min{T 1) T 2 }, 
with Ti ~ (YiUi) 1 / 2 exp(— 2r,/7r|ej|) — are obtained from 



the non- interacting (free-fermion) ones (Table [TTJ) via fu 



TABLE II: Gluing conditions for free fermions 



Qi 


Q2 


3 


Ising 











I 


1 





1/2 


a 





1 


1/2 


a 








1 


e 


1 


1 


1 


I 


1 


1 





e 



sion with the Ising field a [l6| . This means that an Ising 
conformal tower gets replaced according to the operator 
product expansion (OPE) of a with the primary field 
corresponding to that tower, i.e. a x a = I + e and 
exa = Ixa = a. Since we have a free bosonic theory, 
the energy spectrum can then be read off immediately as 



E — Eq 



7T 



(Qi) 2 + (Q 2 ) 2 +. ? (j + l)+4A 



(6) 



with the quantum numbers Qi,Q2,j, and A constrained 
by the gluing conditions obtained from the ones in Table 
|n]by fusion with a (see Table Hv|) . 

The BCFT approach for the TIKM - as reviewed 
above — can easily be adapted to the TIAM, Eq. 
To take into account the fact that, unlike a Kondo im- 
purity, an Anderson impurity (still in the n% = 1 limit) 
carries charge, we redefine the charges as Qi — > Qi — 1. 
To keep the correct physics, the gluing conditions must 
be modified accordingly, by exchanging the Qi = 1 and 
Qi = conformal towers (corresponding to a 7r/2 phase- 
shift of the electrons). These replacements are formally 
equivalent to fusion with the ji = 1/2 conformal towers 
or with the Ising field e in the coset picture. As pointed 
out in Ref. Q, the critical point is symmetric with re- 
spect to 7r/2 phase-shifts, which is formally reflected in 
the OPE a x e = a. Thus the fusion that gives rise to 
the proper boundary conditions is not affected by our re- 
definition of charge. The operator content at the critical 
point is obtained via double fusion [l7| with a x e = a, as 
in Ref. 0], and one finds five operators with dimension 
A < 2 (Table llllj) that are allowed by conservation of 
charge and total spin. 



In the Kondo limit, where particle-hole symmetry is 
present and when there is parity symmetry (^1 , d\ <-» 
1P2 , &2 ) 1 only the first and last operators are allowed 0] • 
The operator e is the only relevant one and comes with a 
scaling field that measures the distance of the RKKY cou- 
pling K to its critical value. Breaking parity symmetry 
will only allow the irrelevant operator Jf_ 1 ■(/>, which docs 
not drive the system away from criticality. If we leave the 
Kondo limit and allow charges to fluctuate between the 
dots and the leads, i.e., we have n, ^ 1, particle- hole 
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TABLE III: Operator content (A < 2) of the TIKM at crit- 
icality. is the j = 1 spin-boson field. As we will see, the 
same list of operators and corresponding dimensions describes 
the critical TIAM. 



Operator 


A 


t 


1/2 


Jt ~ tilt 


1 




3/2 


JU ■ cj> 


3/2 


L-ie 


3/2 



symmetry is no longer present and the second and third 
operators in Table Mil become allowed. None of these 
operators is relevant (A < 1) and, therefore, criticality 
is not destroyed. The exactly marginal charge currents 
JP extend the TIKM critical point to a critical surface 
parameterized by the impurity valences rij. The values of 
these arc determined by the Friedel-Langreth sum rule, 
which relates the charge rii of a local perturbation to 
the electronic scattering phase shift, Sfl = nrii/2, at the 
Fermi level. The finite-size energy spectrum is 



E = E ° + Tl 



^(Q,-n i ) 2 +j(j + l) + 4A 



(7) 



where the phase shifts are imparted by the exactly 
marginal charge currents J?, and the quantum numbers 
QiiQ2,j, and A are constrained by the same gluing con- 
ditions as in the case of the TIKM (Table HV]). 



TABLE IV: Gluing conditions at criticality 



Qi 


Q 2 


3 


Ising 











a 


1 





1/2 


I 





1 


1/2 


I 


1 





1/2 


e 





1 


1/2 


e 








1 


a 


1 


1 


1 


a 


1 


1 





a 



The impurity valences rij contribute to the spectrum 



as E{m) = {■K/Ai)Y li= i 



ZniQi)- The quadratic 



term is identical for all states, regardless of their quan- 
tum numbers and just offsets the ground state energy Eq. 
The linear part shifts the energy of states with different 
charges relative to each other. This shift is of the same 
order as the energy gaps between the different states and 
changes the spectrum qualitatively. The lowest energy 
excitation in the TIKM limit for instance becomes de- 
generate with the second lowest at n\i = 3/4 and has 
higher energy for even lower values. 



III. THERMODYNAMICS AND TRANSPORT 

In contrast to the finite-size energy spectrum, the scal- 
ing dimensions of the operators which determine thermo- 
dynamics and transport do not depend on the impurity 
valences. This is because the scaling dimensions of these 
operators related to the energy levels of the finite system 
with the same boundary condition on either side of the 
strip, and are, hence, obtained via double fusion (l7j . As 
explained in Ref. [la ], the replacement Qi — > Qi — rii 
at the lower boundary corresponds to the replacement 
Qi — * Qi + n-i at the upper boundary, thus the contribu- 
tions from both boundaries cancel each other and there is 
no n^-dependence. The list of operators and correspond- 
ing dimensions as given in Table Mil is therefore valid for 
the TIAM with arbitrary values of the impurity valences, 
< m <l(i= 1,2). 

The low-temperature impurity thermodynamics, as 
well as transport properties due to the presence of the 
impurity, are determined by the dimension of the leading 
irrelevant operators, which for the present case are listed 
as the last three operators in Table Hill The first of these, 
Jfe, is a novel composite operator that appears in this 
model. In the parity symmetric case it governs both the 
impurity thermodynamics and finite temperature impu- 
rity corrections to electron-electron correlation functions. 
The next operator in the list, Jt\ ■ </>, is only allowed if 
parity (ipi,d\ <-> ^2,^2) is broken. In this case it pro- 
duces the same scaling as known from the overscreened 
two-channel Kondo model (TCKM) The last op- 

erator in the table, L_ie, is a Virasoro first descendant 
and does not give any contribution to finite-temperature 
properties, contrary to what is claimed in Ref. [8| (in con- 
text of the TIKM) [l9j]. The presence of the composite 
Jfe operator, and — in the case of broken parity symme- 
try — the operator J s _ 1 ■ <fi, produces non-analytic scaling 
of observables, signaling a non-Fermi liquid ground state. 
The scenario is similar to that of the TCKM, which is well 
covered in the literature. An in-depth discussion in the 
same language as we use here can be found in Ref. [l2l |. 
In particular, the first correction to the zero-temperature 
conductance Go across each dot, generated by Jf e, scales 
as 



5G(T) - y/T/TK, 
and the contribution to the specific heat is 
SC(T) - T\n(T K /T). 



(8) 



(9) 



The impurity susceptibility is not affected by Jfe, since 
these operators do not contain any of the SU(2)2 spin- 
boson fields that could couple to a magnetic field. How- 
ever, when parity symmetry is broken, e.g. by tuning the 
tunneling rates such that V\ 7^ V2, the operator Jt\ ■ <fi 
comes into play, giving rise to a logarithmically divergent 
term in the susceptibility, 



^(T)~ln(T/T^), 



(10) 
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thus replacing the "screened" behavior, x(T) = const. + 
0(T), in the parity-symmetric case. 

Let us conclude this section by pointing out that while 
the impurity valences rij do not affect the scaling to any 
order, they are expected to appear in the non-universal 
amplitudes that multiply the expressions in Eqs. (JSJ - 
(|TU|) . To pinpoint the exact dependence of these ampli- 
tudes on rij, however, requires methods beyond a BCFT 
approach [2C| . 

IV. BUILDUP OF THE CRITICAL GROUND 
STATE: ASYMPTOTIC SCALING 

We now shift focus, and consider a situation where 
one or several of the tunnel junctions in the device are 
"pinched off" (by a change of the gate voltages that 
control the corresponding junctions), implying that the 
electrons in the two leads become dynamically decou- 
pled. At a later time — after the severed system has 
relaxed to equilibrium — the device is reconnected by a 
sudden "switch on" of the same tunnel junctions (again 
by a change of gate voltages). This poses a question 
that is a variant of the Fermi edge singularity problem: 
"How much time does the system need for rebuilding the 
non-Fermi liquid critical state after the sudden pertur- 
bation?" . To answer this question we consider a setup 
initially described by two decoupled Fermi liquids, the 
exact nature of which depends on which tunnel junctions 
are considered. Then, at some time t = to we suddenly 
reconnect the system by turning on these tunnel junc- 
tions. If the parameters of the setup are properly tuned, 
the system will evolve into the quantum critical phase 
associated with the TIAM non-Fermi liquid ground state 
(Fig. ®. 



Lead 1 Lead 1 Lead 1 




Lead 2 Lead 2 Lead 2 



FIG. 2: (Color online) One of the physical processes that we 
consider in this section (case Ci in Table V). In the first step 
there is a SIAM and a decoupled channel of free electrons. 
In the second step, the coupling Vi is suddenly turned on, 
allowing charge to fluctuate between the upper lead and the 
upper dot. In the last step the system has relaxed to its 
critical ground state due to the RKKY interaction between 
the two dots. 

We now systematically consider the different scenar- 
ios that can emerge, corresponding to turning off (and 
on) the various (combinations of) couplings (see Tabk fV|) . 
Before we proceed with details, let us briefly discuss the 
various possible processes. Process "A" is maybe the 
most obvious one: If both couplings between the leads 
and the dots are turned off, the initial state of the con- 



TABLE V: List of "buildup" processes. The boldface crosses 
denote couplings that are first turned off, and then — after 
the severed system has equilibrated — are suddenly turned 
on again. Taking 1 <-> 2 completes the list to all possible 
combinations. The crosses in normal font mark couplings that 
may or may not be turned off and on without changing the 
character of the respective process. 







v 2 


vl 


vl 


initial lead configurations 


A 


X 


X 


X 


X 


1. free electrons 

2. free electrons 


B 






X 


X 


1. SIAM 

2. SIAM 


d 


X 




X 


X 


1. free electrons 

2. SIAM 


c 2 


X 








1. free electrons 

2. free electrons 



duction electrons is that of a free electron gas. Note 
that there are actually three different variations of this 
process, namely with both (decoupled) dots empty, with 
one of them empty and one of them filled, or with both 
of them filled. Process "B" , i.e. turning off one of the 
dot-reservoir couplings, say V~a,i, turns the system into 
a SIAM for lead 1 and its attached dot, while the rest of 
the system gets described by an anisotropic two-channel 
SIAM (with the two channels provided by lead 2 and the 
auxiliary reservoir, respectively). Since the latter model 
renormalizes to a single-channel SIAM [22j the conduc- 
tion electrons are effectively described by two decoupled 
SIAMs, just as if both Va.i and Va.2 were turned off. 
The last two processes in Tabic V, Ci and C2 — both 
corresponding to turning off the coupling V\ — are a bit 
more involved and depend on whether the dot 1 carries 
charge or not after decoupling from lead 1. Certainly, 
the conduction electrons on lead 1 are described by a free 
electron-gas in either case. If the dot 1 carries no charge, 
then lead 2 is described by an anisotropic two-channel 
SIAM which renormalizes to a single-channel SIAM, as 
before. On the other hand, if the dot is charged, it 
will form a singlet with dot 2 since, by the constraint 
K ;§> Ta, the RKKY interaction dominates the Kondo 
interaction. As a consequence, the second lead decouples, 
turning it into a free electron gas. 

In the language of BCFT the interactions with the im- 
purities only appear via the boundary condition, and 
hence the physical process of changing couplings cor- 
responds to changing the boundary conditions. This 
can formally be described by an appropriate boundary- 
changing operator acting on the system. To make this a 
bit clearer, we introduce a notation as depicted in Fig. [3] 
for each of the two possible channels that may enter into 
the initial configuration (the severed system at equilib- 
rium): free electrons [isolated lead] and SIAM [lead with 
an attached dot] . Since the severed system consists of two 
leads (with, or without, an attached dot) both must be 
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specified in order to determine the initial configuration. 
We here note that for a decoupled system, the fusion and 
charge redefinition that take the free electron gas onto a 
S1AM can be performed successively. The trivial fusion 
with the identity operator I (cf. FIG. 3) reflects the fact 
that the low-energy physics of the SIAM is that of a local 
Fermi liquid, hence the gluing conditions on the charge- 
and spin sectors are the same as for free electrons. 



Free electrons 
SIAM 



FIG. 3: (Color online) Notation for the single Fermi-liquid 
leads (i — 1, 2) and corresponding fusion/redefinitions. Here 
■Krrii/2 is defined as the phase shift that the electrons expe- 
rience at x = 0, such that mi equals the net charge on the 
dot (d\di) when it is coupled to the leads (and m,i = when 
the coupling is turned off, independent of {d\di)). The sign 
in front of nn depends on whether the dot is attached to the 
lower (+) or upper boundary (-) of the finite strip. 

The charge redefinition Qi — * QiT mi follows as before 
from the Friedel-Langreth sum rule. A remark is in order 
concerning rm, which is the net charge on dot i when it 
is coupled to lead i but decoupled from the other dot. It 
does not have to be equal to rii, which is the net charge 
of the same dot at the critical point of the TIAM (with 
the two dots being RKKY-coupled) . In fact, we expect 
that rii < rrii {ni,mi 6 [0, 1]). Since only the singly oc- 
cupied subspaces participate in the RKKY interactions, 
the effective RKKY coupling increases with increasing 
rii. The RKKY coupling is antifcrromagnctic (K > 0) 
and the spin-spin expectation value at the critical point 
is negative [23| , thus, larger values of are energetically 
favorable. 



Two-impurity 
Anderson Model 



x<r, Qi 



>1,2 T n l,2 



FIG. 4: (Color online) Notation for coupled leads and corre- 
sponding fusion/charge-redefinitions. 

We represent a decoupled system by drawing the two 
leads parallel to each other, using the notation from FIG. 
3. The two leads can be connected only by the RKKY 
interaction, which we denote by a line connecting the 
impurities (Fig. [4j. Using this notation, we can represent 
the process of turning on the couplings (lead-dot tunnel 
junctions and/or the RKKY interaction) at r = r = it 
as shown in Fig. [5J where the initial configuration is 
chosen here as two SIAMs (entry B in Table V). The 
quantity that measures how far the system has evolved 
from its initial state after a time At, is given by the 
Green's function of the (yet undetermined) boundary- 
changing operator O: 



G(t) = (0,I|O(r + Ar)Ot( r0 )|0,I) ~ (At) 



(11) 



where x is the scaling dimension of O and |0,I) is the 
initial ground state of the severed system. 



r = r 

FIG. 5: (Color online) A process where the initial state con- 
sists of two decoupled SIAMs, where the boundary condition 
is changed at tq to that of the TIAM. 



The boundary-changing operator is assumed to be a 
primary field so that the asymptotic states are ground 
states of their respective configurations. In this case the 
Green's function is slaved to the given asymptotic form 
by conformal invariance. To relate x to the energy spec- 
trum, we map the half plane onto the infinite strip via 
the conformal transformation u + w = (1/tt) ln(r + ix) 
(Fig. 0, following the prescription laid out in Ref. [2l| . 

t t? e 1 



r = 



1= 



FIG. 6: (Color online) Mapping to the strip geometry results 
in a theory with boundary condition I on one side and bound- 
ary condition F on the other. 

Under this transformation the propagator takes the form 



(21 7T 

(0, H|0(w o + Au)O f {u a )\0, II) = I — sinh — Au 

A«»( ™ -ZxAv 
I 



(12) 



where uq = (d/ir) Ln(ro). Here we denote the ground 
state of the system on the strip geometry (same bound- 
ary condition on both sides) by |0, II). Alternatively, the 
propagator can also be evaluated in the strip geometry 
by inserting a complete set of states. Note that since O 
changes the boundary conditions, the eigenstates of the 
system after applying O are not the same as the ones 
before (see Fig. [6]). We label states with boundary con- 
dition I on the top and boundary condition F on the 
bottom (see right hand side of Fig. [6]) by |n, IF). These 
are the eigenstates after applying O and a complete basis 
of them should therefore be inserted, 



<O,n|0(u o + Au)O t (iio)|O,II) 

= 5^(0, ll\e HAu O(u Q )e- HAu \n, IF) (n, H^MIO, II) 



:^|(0,II|O( W0 )|n,IF)| 2 e A "«-^ F ). 



(13) 
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In the limit of large Au, only the lowest-energy state 
no, IF) contributes and by comparison with Eq. (|12| . we 
get: 



t 



(14) 



So far, no particular boundary-changing operator 
has been specified. Since the scaling dimension of a 
boundary-changing operator depends only on the energy 
of the corresponding state on the strip, with one type of 
boundary condition at the bottom and another type of 
boundary condition on the top, the possible boundary- 
changing operators are determined by the gluing condi- 
tions in this geometry. Using the notation introduced 
above (cf. Figs. [3] and 2]), it is a straightforward task to 
find the correct operators and their scaling dimensions. 
The energy is measured with respect to the ground state 
of the system with boundary condition II, which cor- 
responds to two decoupled SIAMs on each side of the 
strip. According to Fig. [3] the gluing conditions are the 
ones for free fcrmions ( Table IIT|) . Moreover, the redefi- 
nitions of the charge on either side of the strip cancel: 
Qi Qi + frii — rrii = Qi . Thus Eq 1 is indeed equal to 
the ground-state energy of the free fermion theory Eq. 

Since the SIAM connects continuously to the free elec- 
tron gas via the charge redefinitions Qi^rrii it is sufficient 
to consider two decoupled SIAMs as the initial state and 
get the other possibilities by setting the appropriate m, 
to zero and keeping track of the correct quantum num- 
bers Qi, Qi. For the configuration labeled by IF (Fig. [51 
right hand side), the boundary condition is that of two 
decoupled SIAMs at the top and the one of the TIAM 
at the bottom of the strip. By comparison with Figs. [3] 
and @] this translates into fusion with a and the charge 
redefinitions Q\ — > Q\ — n\ + m%, Qi — > Qi — n<i + mi. 
Defining hi = rii — m,; the energy for states with these 
boundary condition is 



E = E ° + Tl 



^(Q,-^) 2 +X? + l) + 4A 



(15) 



Combining this with Eq. (fFlj) immediately yields the 
exponents for the various processes (see Table IVlf . 



TABLE VI: Exponents governing the buildup of the quantum 
critical ground state for the various processes. To distinguish 
between the different possibilities for each process (see dis- 
cussion following Table [V| , we have specified the initial net 
charge on the dots ({d\di), {d^d^)) for each process. 



Process 


Xi 


A(0,0), B(mi,m 2 ), Ci(0,m 2 ) 


i + \h\ + \h\ 


A(1,0), Ci(l,m 2 ) 


- \h x + \h\ + \h\ 


A(l,l), G a (l,l) 





The Green's function of the three different bound- 
ary changing operators 

Gfj(t) = (O|O<(t)0t(o)|O) ~ i = 1,2,3 

measures the overlap between the initially perturbed 
state and the state to which the system has evolved at 
time t, with Xi listed in Tab (VI). It is interesting to 
compare these exponents with the ones for two decou- 
pled SIAMs, i.e., the process depicted in Fig. [Jj where 
the fiij of the two SIAMs take the role of the hi for the 
TIAM. The analogy is most obvious in the cases A and 
Ci when hi = n,;, i.e., when the initial state consists of 
two channels of free electrons, (as in Fig. [7j). 



r = r 

FIG. 7: (Color online) Two decoupled channels of electrons 
are turned into two decoupled SIAMs at to. 

Since the two SIAMs are completely independent, the ex- 
ponents for two such models [24j | can simply be added and 
we find that our exponents in Table VI (upon identifying 
mi and hi) differ from those only by constants of ±1/16, 
where 1/16 is the scaling dimension of the Ising field a. 
The fact that the valences guide the respective processes 
in this very similar way is notable, considering that the 
TIAM critical ground state is of non-Fermi liquid type. 



V. SUMMARY 

In conclusion, we have proposed a realization of the 
two-impurity Anderson model in a double-quantum dot 
device with separated leads. We have presented an exact 
BCFT description of the system at quantum criticality 
and used it to calculate several physical properties. We 
have found a new composite operator that governs the 
low-temperature behavior, and explained how non-Fermi 
liquid scaling behavior which is not present in the original 
TIKM Q emerges as a consequence of certain (irrelevant) 
symmetry breakings. Whereas the energy spectrum ex- 
plicitly and qualitatively depends on the charge valences 
of the dots, thermodynamics and transport properties are 
only affected in non-universal amplitudes. Instead, the 
charge valences enter in the dynamical exponents that 
set the time scale for buildup of the non-Fermi liquid 
quantum critical state. 
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